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ABSTRACT 

10 ' Context. Long-lived, large-scale magnetic field configurations exist in upper main sequence, white dwarf, and neutron 

stars. Externally, these fields have a strong dipolar component, while their internal structure and evolution are uncertain 
but highly relevant to several problems in stellar and high-energy astrophysics. 

Aims. We discuss the main properties expected for the stable magnetic configurations in these stars from physical 
arguments and the ways these properties may determine the modes of decay of these configurations. 
Methods. We explain and emphasize the likely importance of the non-barotropic, stable stratification of matter in all 
^ , these stars (due to entropy gradients in main-sequence envelopes and white dwarfs, due to composition gradients in 

5_j . neutron stars). We first illustrate it in a toy model involving a single, azimuthal magnetic flux tube. We then discuss the 

effect of stable stratification or its absence on more general configurations, such as axisymmetric equilibria involving 
' poloidal and toroidal field components. We argue that the main mode of decay for these configurations are processes 

that lift the constraints set by stable stratification, such as heat diffusion in main-sequence envelopes and white dwarfs, 
and beta decays or particle diffusion in neutron stars. We estimate the time scales for these processes, as well as their 
CO ' interplay with the cooling processes in the case of neutron stars. 

^ ' Results. Stable magneto-hydrostatic equilibria appear to exist in stars whenever the matter in their interior is stably 

stratified (not barotropic). These equilibria are not force- free and not required to satisfy the Grad-Shafranov equation, 
. but they do involve both toroidal and poloidal field components. In main sequence stars with radiative envelopes and 

CO ' in white dwarfs, heat diffusion is not fast enough to make these equilibria evolve over the stellar lifetime. In neutron 

' stars, a strong enough field might decay by overcoming the compositional stratification through beta decays (at the 

highest field strengths) or through ambipolar diffusion (for somewhat weaker fields). These processes convert magnetic 
energy to thermal energy, and they occur at significant rates only once the latter is less than the former; therefore, they 
substantially delay the cooling of the neutron star, while slowly decreasing its magnetic energy. 
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1. Introduction $max = ttR'^B^t.^^ ^ 10^^-^ G cm^. This coincidence has 
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JlJ . . 1 • 1 r 1 often been interpreted as an argument for flux freezing 

. .Upper mam sequence stars, white dwarfs and neutron stars during stella r evolution ( Rudcrman. 1972; RciscncW 

appear to have long-liyed magnetic fields. These fie ds are [2003I ; iFerrario fc Wickramasingha l2005allU I200i) . 

orgamzed on large scales, m the sense that the dipole field although its feasibilit y has been called into question 

com ponents (and pe rhaps some other, low-order multipoles; (|ThomDson fc Duncanl . \l99^ Is^ . [200I . Of course, a 

e g. i Bagnulo et alj^999 lOOa) are not much weaker than ^ ^^^^^.^^ ^^.^^^^ magnetic flux might be ejected 

the rms surface fleld, unlike the highly chaotic field of the ^-^^ g^^^^^ envelope. On the other hand, substan- 



S 



un. 



1 / r 1- 1 \ -nil tial field amplification through differential rotation, con- 

The highest detected (surface dipole) magnetic field ^^^^^^^ ^^^.^^^ instabilities could plausibly occm 

3 X 10 G m chemicallly peculiar A and B stars (Ap/Bp Connected to the similar fluxes is that these stars also 

stars , R ^ 3 Rq; | Mathys et al. | ll997t i Bagnulo et al. | ^ave similar ratios of fluid pressure {P ~ GM'^/R^ where 

19991) ■ Bm.^ ~ 10 ^ G in white dwarfs (i? ~ lO^'km; G is the gravitational constant and M is the mass of the 



Schmidt etlH l2003f ). and Bmax ^ 10^^ G in "magne- star) to magnetic pressure (S^/Stt), 
tars", a sub class of strongly m agnetized neutron star s 

{R ~ 10 km; iKouveliotou et al.l [r998; Wo ods et al.lll999tl . „ _ SttP Sn^GA'P ind f V f ^ 



in all cases yielding very similar total magnetic fluxes, ^ B'^ ^ $2 '^-'3x10 ^^^^J ' ' ^^"^ 
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a large number, even for the most highly magnetized ob- 
jects, implying that the magnetic field causes only very mi- 
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nor perturbations to their hydrostatic equilibrium struc- 
ture. 

Another similarity among these stars is that much or 
all of their structure is stably stratified, i. e. stable to con- 
vection. The radiative envelopes of upper main sequence 
stars, as well as the whole interior of white dwarfs, are sta- 
bilized by the radially increasing specific entropy s, while 
in neutron stars the same effect is caused by a r adially 
varying mix of different particle species (Pcthick, Il992t 
iReisenegger fc Goldreichl Il992t iReiseneggCT , ,2001a) , which 
in the outer core reduces to a radially varying proton and 
electron fraction, Y = np/(nn + Up) = ne/(n„ -I- rip), where 
rii stands for the number densities of neutrons {i = n), 
protons (i = p) and electrons (i = e). 

The structure of the magnetic field inside these stars is 
not known, although it is highly relevant to their evolution: 

1) It affects the radial trans port of angular mo mentum and 
chemical elements (e. g. iHeger et al.|[2005l ): 

2) it is plausibly the dominant source of energy for 
both the outbursts and the persistent emission of soft 
gamma repeaters (SGRs) and anomalous X-ray pulsars 
(AXPs), for this reason collectively called "magnetars" 
(^Thomp son fc Duncanl . Il993l . 119961) : 

3) it is likely to play an important role in the frequency 
spectrum o f quasi-perio dic oscillations observed after 
SGR flares (Levin|, 12007D; 

4) it leads to slight deformations of neutron stars that 
could give rise to precession of pulsars ( Wasserman, 
2003i) and to the emission of gravitational waves ((Cutler , 



Magnetohydrodynamic (MHD) simulations of stably 
stratified stars with random initial magnetic field configura- 
tions have shown them to evolve on Alfven-likc time scales 
into long-lived structures whose further evolution and de- 
cay appears to be controlled by dissipa tive processes (in the 
simulations, Ohmic diffusion; Braith waite fc Spruit] |2004 
12006': 'Braithwait e fc Nordlundl[2006il . Often, but not always 
(^Braithwaite, 2003), these are roughly axisymmetric combi- 
nations of linked poloidal and toroidal components, whose 
external appearance is essentially dipolar. It appears plau- 
sible that these configurations approximate the true mag- 
netic field structures in upper main sequence stars, white 
dwarfs, and neutron stars. 

In § 2, we present arguments to the effect that 
the stable stratification of the stellar matter should be 
an essential ingredient to these equilibria. This means 
that, contrary to assumptions in the recent literature 
(e. g. Perez- Azorin et al. 20061 : iBroderick fc Naravan 2008i : 



iMastrano fc Melatosll2008l ). these are definitely not force- 
free fields. In fact, it is shown in Appendix A that 
there are no true force-free equilibria in stars, while 
those proposed in the literature actually have singu- 
lar magnetic forces on the stellar surface. Moreover, 
the fluid cannot be treated as barotropic, therefore the 
field components arc not required to satisfy the Grad- 
Shafranov equation (Mcstcl, 1956), contrary to the pop- 
ular belief (Tomimura & Eriguchi, 2005; Yo shida et al.l . 
20061 [Haskell et al... ,2008; .Akgiin fc WassernlMl 120081: 



Kiuchi fc Kotakel . l2008ir In fact, the range of available 



equilibria becomes much wider in a stably stratified, non- 
barotropic fluid. The constraints imposed by the stability 
of these equilibria are far from obvious, but we argue that 
there are probably no equilibria in barotropic stars, while 



it is likely that there are equilibria with linked toroidal and 
poloidal flelds in stably stratifled stars. 

Of course, the specific entropy s and the proton frac- 
tion Y are not perfectly conserved quantities within each 
fluid element, but can be changed by dissipative pro- 
cesses, dis cussed i n 3 : In the case of s, through heat 
diffusion (jParkerl |1974[ ). in the case of Y, by (direct 
or inverse) beta decays or by ambipolar diffusion (mo- 
tion of charged relative to neutral particles; Pethick 1993; 
Goldreich fc Reiseneggeijri992l : lThompson fc Duncanlll996l : 
Hovos et al.l 20081) . Thus, the condition of stable strat- 
ification, and with it the hypothetically associated sta- 
ble magnetic equilibrium configuration, although excel- 
lent approximations on short (Alfven-like) time scales, are 
eroded on the time scales of the dissipative processes men- 
1 :ioned above, leading to the d ecay of these structures 
([Goldreich fc Reiseneggeii Il992f). and perhaps to a sud- 
den loss of stability (lBraith\vaite fc Spruitl . |2004 [20061 : 



[Braithwaite fc Nordlundl . [20061) . In" main-sequence stars, 
white dwarfs, and weakly magnetized neutron stars, these 
appear to be too long to act on the stellar life time, but 
in strongly magnetized neutron stars their time scales be- 
come shorter, so they might plausibly drive magnetic field 
decay, lea ding to internal heating and to the magnetar phe - 
nomenon ([Thompson fc Duncad . [I996t [Arras et al.l [2004[ ). 
In principle, the Hall drift might also play a role 



in the field evolu t ion in neutron stars (iJonesl. 



1992 



Urpin fc Shalvbkovl 1991*: ^G oldreich fc Reiseneggeii . _ 
Reisenegger et al.l [2005. 2007[: [Pons fc Gepperfj . [2007[ ) 



1988 



However, its time scale in neutron star cores tends to be 
somewhat lo nger than those of the other processes con- 
sidered here ([Goldreich fc Reiseneggeilll992[ ). Moreover, in 
an axially symmetric, equilibrium magnetic field configu- 
ration, the effect o f the H all drif t is exactly cancelled by 
bulk fiuid motions ([Reisenegger fc Thompsod . [2009') , so we 
do not take it into account. For simplicity, we also refrain 
from discussing the role of the solid crust of the neutron 
star, as well as the effects of superconducti vity and superflu- 
idity, which alter the magne t ic str esses (E asson fc PethicS . 
[igTTt [Akgiin fc WassermM . [20081) as well as the dissipa- 
tive processes. We also ignore the process of initial set-up 
of the magnetic equilibrium, which might be a highly dy- 
namical process involving differential rotation and possibly 
a dynamo (Thompson & Duncan, 1993; Spruit, 2002), but 
concentrate on the properties imposed by the equilibrium 
and stability conditions and on the long-term evolution of 
the field. 

A concise summary of our conclusions is given in § 4. 
Parts of this discussion h ave already been given elsewhere 
([Reiseneggeii [20071 [20081 ) . 



2. Magnetic equilibria and stable stratification 

2.1. Force balance 

In a conducting, fluid star, a general MHD equilibrium is 
set by the condition that the net force on the fluid vanishes 
everywhere, i. e. 



f B + fF = 0, 

where 

fs = ij X B 

c 



(2) 



(3) 
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is the magnetic (Lorentz) force per unit volume, written in 
terms of the magnetic field, B, and its associated current 
density, j = (c/47r)V x B, and 



ip = -VP - pS/ip 



(4) 



is the fluid force, which depends on its pressure, P, density 
p, and gravitational potential, ip. 

In all the stars of interest, the fluid is non-barotropic, 
i. e. the density is not a function of pressure only, but de- 
pends on an additional, non-trivial variable X, which is 
conserved on dynamical (sound or Alfven wave crossing) 
time scales: specific entropy {X — s) in the case of up- 
per main sequence stars and white dwarfs, and the fraction 
of protons {X = Y) or other minor constituent particles 
required by beta equihbrium in the case of ne utron stars 
(|Pethickl [l99llReisenegger fc Goldreichl . II992D . 

As shown in equation ([T]), the fluid pressure is much 
higher than the magnetic pressure, so we take the point 
of view that the magnetic forces create only a slight per- 
turbation to the background hydrostatic equilibrium state 
(ip = 0) the star would have in their absence. In addition, 
we invoke the standard Cowling approximation of neglect- 
ing perturbations to the gravitational potential, also used 
in the simulations of Braithwaite and collaborators. 

We do not assume that the unperturbed star is spher- 
ically symmetric, so our arguments can also be applied to 
stars that are uniformly rotating, in which case tjj has to 
be interpreted as the effective potential, also including cen- 
trifugal effects. However, we ignore the effects of merid- 
ional circulation. The time scale for this process, due to 
the interaction of stellar rotation and internal heat flow, is 
tchc ^ {^K /^)^tKH, where is the stellar rotation rate, 
Hk is its maximum (Keplerian or "break-up") value, and 
^KH is the (Kelvin-Helmholtz) time scale required to ra- 
diate away the thermal energy content of the star. For 
main sequence stars, ^kh is substantially shorter than their 
main-sequence life time, so meridional circulation can mod- 
ify the magnetic field s tructure of suffici e ntly fas t rota - 
tors, P = 2TT/n < 5d (iMossl . 11984 llQQfll : iMestel Il999t) . 
to which our analysis will therefore not apply. For white 
dwarfs, ^kh is their cooling time, i. e. essentially their age, 
so meridional c irculation is unimpo rtant unless they rotate 
near break-up (jTassoul fc Tassoul . Il983f ). In none of these 
cases, the magnetic field is expected to be strong enough 
to have a substantial influence on the pattern or time scale 
of meridional circulation. For neutron stars, since the main 
source of stratification is not entropy b ut chemical compo- 
sition ([Reisenegger fc Goldreichl . Il992f ) , meridional circula- 
tion will not occur at all. 

Thus, we write the fluid force in terms of the Eulerian 
perturbations (changes at fixed spatial positions) of density 
and pressure, Sp and SP, respectively, as 



(5) 



The perturbations can be viewed as being produced by 
a displacement field ^, which allows us to introduce 
Lagrangian perturbations (comparing the variables in the 
same fluid element as it gets displaced) Ap, AP, formally 
related to the Eulerian perturbations by 



(6) 



where the gradient operator acts on the corresponding un- 
perturbed "background" quantity, pf, or Pi,. 



If the displacement is fast enough, it can be taken to 
conserve specific entropy and chemical composition, so the 
Lagrangian perturbations are related by the adiabatic index 



7p 



AlnP 
Aln/9 



ainP 
d\np 



(7) 



X 



which is generally different from the analogous quantity 
characterizing the hydrostatic equilibrium profile of the 
star. 



dlnPh 



lb 



dlnpb \ 9 In p 



ainP 



ainP\ dlnX 



dlnX J dlnph 



(8) 



Then, also using the condition of mass conservation, 
Ap — —pb^ ■ ^, ec^uation ^ can be manipulated into the 
form 



(9) 



where the first term would be present as well in a 
barotropic, homogeneous fluid, whereas the second ac- 
counts for buoyancy effects. The latter is stabilizing if 
7p > 76 s-nd destabilizing in the opposite case. 

In upper main sequence stars, the fluid is a clas- 
sical, monatomic ideal gas with 7p = 5/3, with 
their radiative envelope s well described by 7^ ~ 4/3 
(|MacGregor fc Cassinellil . I2003D . so 7p/7f, - 1 « 1/4. In 
white dwarfs, the electrons are highly degenerate {kT <C 
Epe, where k is Boltzmann's constant, T = x lO^K is 
the interior temperature, and Epe is the electron Fermi 
energy, not including the relativistic rest-mass term, mc^) 
and dominate the pressure, but the entropy is contained 
in the ions, so 7p/7h - 1 ~ kT/ZEp - T7/5OO. Finally, 
in the case of neutron stars, entropy becomes negligible 
a few seconds after their birth, but they remain stably 
stratified due to the den sity-dependent proton fr a ction 
7D/ 7b - 1 ^ Y few % (iReisenegger fc Goldreichl 119921: 
iLaJl . 119941: iReiseneggeii ^2001a^ . 

Eq. ^ shows that, in a stably stratified fluid (with 
7p > 76); the fluid force has two parts that are deter- 
mined by two independent, scalar functions, e. g. 6P and 
Ap, which give the fluid a greater freedom to balance the 
magnetic force than it would have in the barotropic case 
ilp — lb)- It should be noted that, if the buoyancy term 
in eq. ^ is crucial to balance a particular field configura- 
tion of characteristic length scale comparable to the stellar 
radius, the characteristic field strength is constrained by 



< 



SttP 

lb 



^-1) lApV^I 
lb J 

1 ) pIvv^I « J - 



1 IVPI 



(10) 



Thus, its maximum value is not set by the condition /? > 1 
(with /3 defined in eq. [T]) , but rather by the more restrictive 



Pilvllb - 1) > 1, 



(11) 



yielding maximum allowed field strengths of ~ 10* G for 

Ap/Bp stars, ^ lO^^T^^^ G for white dwarfs, and a few 
times 10^^ G for neutron stars, all of which still substan- 
tially exceed the observationally inferred flelds in these 
stars. 
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2.2. Hierarchy of equilibria and variational principles 

Given this physical background, we now explore the more 
mathematical issue of how magnetic equilibria obtained un- 
der progressively more stringent (and, for our purposes, 
more realistic) constraints can be represented by con- 
strained stationary points of the magnetic energy. We note 
that, contrary to the previous section, where "perturba- 
tions" were deviations from the non-magnetic, hydrostatic 
equilibrium caused by the magnetic field, here perturba- 
tions are taken with respect to successive, magnetic equi- 
libria. 



2.2.1. Field-free: 

Consider the total magnetic energy within a fixed volume, 
Ub = J^B'^dV/iSir). The only way to obtain SUb = 
under a weak, but otherwise arbitrary magnetic field vari- 
ation is to have B = everywhere. This is the absolute 
minimum of the magnetic energy, and it is eventually ob- 
tained in a star placed in vacuum, without external fields, 
and in which a sufficiently effective dissipation mechanism 
(such as resistive diffusion) is active. 



2.2.2. Current-free: 

Of course, (5B is not fully arbitrary, but must be diver- 
genceless, so we now consider the slightly more restricted 
case (5B = V X i5A. Now, 



SUb = — / B-V xSAdV 
47r Jv 



1 

An 



I V • ((5A X B) dV + / (5A • V X B dV 



(.12) 



In the last result, the first integral can be made to van- 
ish by imposing appropriate boundary conditions on the 
surface, for example that the magnetic flux through any 
surface element is fixed (n • (5B = 0, where n is the outward 
unit normal), which implies that fi x i5A — 0. The second 
term can only be made to vanish for arbitrary (5A if the 
electric current density vanishes everywhere in the volume, 
j = cV X B/ (47r) = 0. In order to have a finite magnetic field 
in a volume containing no currents, there must be source 
currents in some neighboring volume. Therefore, a finite 
stellar magnetic field cannot be current-free everywhere. 

2.2.3. Force-free: 

In stars, dissipative processes such as the resistive diffusion 
of magnetic flux are slow (see § [3] for details) and therefore 
often negligible. In this context, the only possible perturba- 
tions of the magnetic field arc those which can be produced 
by plasma displacements, 5 A = ^ x B, with ^ an arbitrary 
vector field. In this case, neglecting the divergence term of 
eq. (fT2)) . 5Ub = — /y € • (j ^ B/c) dV, so the stationarity 
of the magnetic energy implies the vanishing of the Lorentz 
force, j X B/c = 0. This case is relevant in very diffuse 
plasmas such as those in stellar magnetospheres, where the 
conductivity is high, but gas pressures and densities are low 
(/3 <^ 1), so the dynamics is dominated by the magnetic 
field. However, as shown in Appendix [XI it is not possible 
to have magnetic field structures that are force-free every- 



where in a star, unless it is confined by (unrealistic) surface 
forces. 



2.2.4. Force balance: 

As we saw in §[231 the plasma inside stars has /3 ^ 1, so the 
fluid forces, due to pressure and gravity, can by no means 
be neglected. The total energy perturbation (with respect 
to some hydrostatic or hydromagnetic equilibrium state) 
caused by an arbitrary fluid displacement field ^ can be cal- 
culated by integrating the work per unit volume — (l/2)f 
done against the forces of Eqs. ^ and (O in order to build 
up this displacement field. In the absence of a magnetic 
field, the result can be written as 



5U 



1 



5P' 



lb 



1 (I • ^PbY 



(13) 



where we ignored a surface term that vanishes for ap- 
propriate boundary conditions. Clearly, the SP^ term is 
always positive, while the other term is positive for a 
stably stratified fiuid (7^ > 7^), zero for a barotropic 
fluid (7p = 7f,), and negative for a convectively unstable 
fluid (7p < 7fc). Clearly, there are no unstable (negative- 
energy) perturbations if 7p > 7b. In a stable equilibrium, 
one must have vanishing Eulerian pressure perturbations, 
SP = -7pPbV • ^ - I • VPb = 0, and, if the fluid is stably 
stratified, also no "vertical" displacements, ^ • VPh = 0. 
All displacements satisfying these conditions produce neu- 
trally stable perturbations. We expect these conditions to 
still be approximately satisfled in equilibria that involve a 
weak magnetic field. 

When a magnetic field is introduced, the total energy 
perturbation becomes more complicated. 



6U 



SB^ 1 

(I • VP)V • ^ - (I • V^)V • (p^)] dV 



j-JBx| + 7pP(V-|)2 



(14) 



([Bernstein et all . [19581 ). and negative-energy pertur bations 
exist for some field configurations, even if 7p > 7f, ([Tavleil 
_1973; Fl owers fc Ruderman , 1977,). However, since these in- 
stabilities are caused by the magnetic field, we expect the 
quantities \6P\ and (if the fluid is stably stratified) || • VP| 
to be small enough to keep the fluid force perturbation 
Sip = — V(5P — SpWip not much larger than the mag- 
netic force perturbation 5{b — (<5j x B + j x 6B)/c. This 
requires |V • (p^)| < p\^\/{P£) for barotropic and stably 
stratified fiuids, and j^^l < only for the latter. Since 

the non-magnetic terms in Eq. p4p are quadratic in these 
quantities, they will be smaller than the magnetic terms, 
therefore it seems reasonable to consider only magnetic en- 
ergy perturbations, subject to the conditions V • (p^) = 
for barotropic and stably stratified fluids, and additionally 
^ • Wtp = for the latter. Note that this bypasses some 
instabilities that originate in particular regions where the 
fluid fo r ces ar e weak, such as near the center of the star 
(lTavler] . [T97l . 

Barotropic fluid: Thus, for a barotropic fluid with high f3, 
we impose the constraint 6P/"fp = Sp = —V • (p^) = 0, 
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therefore ^ = (l/p)V x a, where a is an arbitrary vector 
field. The magnetic energy perturbation becomes 



SUb= I V • 

'V 



(j X B) X a 



pc 



a • V X 



j X B 

pc 



(15) 



Consistent with the condition that 6p — everywhere, we 
require n • ^ = on the surface of the star, implying n x a = 
0, which makes the first integral vanish. The vanishing of 
the second for arbitrary a requires 



vxi'iilSi.o, 

pc 



(16) 



i. e. that the Lorentz force per unit mass be a gra- 
dient (consistent with it having to be balanced by 
the first term in eq. [5]). This is the case considered 
in most expl icit descriptions of n e utron star magnetic 
fields so far (iTomirnura &z Eriguchil. l2005t lYoshida et al 



20061: iHaskell et all. 120081: lAkgiin fc Wassermanl . I200S : 
Kiuchi fc Kotakd . 120081)" 



Stably stratified fluid: Finally, in the strongly non- 
barotropic, stably stratified case (the most realistic, accord- 
ing to our discussion in § I2.ip , vertical motions are strongly 
suppressed, so ^ also has to be tangent to gravitational 
equipotential surfaces, which is equivalent to requiring that 
a = /VV', where / is an arbitrary scalar function. Now de- 
manding that 5Ub — for arbitrary /, we obtain 



V^.Vx(i^)=0, 

a weaker condition than eq. ((T6|) , that is satisfied if 



= V/i + vVip, 



(17) 



(18) 



j X B 

pc 

where p and v are arbitrary scalar functions, consistent 
with eq. 

It is clear that both the limit of a very diffuse plasma 
(/3 <C 1) applicable to § 12.2.31 and the very dense plasma 
{(3 ^ 1) of the present section are idealizations. A more 
rigorous description would minimize the total energy of the 
star, including internal and gravitational energies in ad- 
dition to the magnetic energy, and not impose the addi- 
tional conditions of this section on the displacem ent field 
^. Moreover, the Tayler instability (|Tavleij [19731 see also 
§ 12.41 of the present paper) , although active for arbitrarily 
high values of /3, requires to account explicitly for these 
other contributions to the energy, and consistently to re- 
lax the constraints on ^. For our purposes, the idealized 
description given here appears to be sufficient. 



2.2.5. A note on helicity conservation: 

A variation of the magnetic helicity, H = J^^ A ■ B dV, can 
be written, aside from a surface term (jSpruitj . 120081 ). as 
SH ^ 2 J^TS ■ SA dV, so (for appropriate boundary condi- 
tions) helicity is automatically conserved if 5 A = ^ x B. If 
we initially allow for an arbitrary SA, but the n search for a 
stationary point of Ub at fixed Ti (as do ne bv'Woltierl ll958l 
and more recently bv iBroderick fc Nara van 2008), we ob- 
tain j = aB, where a is a constant Lagrange multiplier. 
This condition is less restrictive than those obtained in 



§ Upland §12X1 but more restrictive than those of 
and 12.2.41 In particular, the force-free condition of 12.2.31 
allows for j = Q;(r)B with B • Va = 0, so a is constant 
on a given field line, but possibly different on different field 
lines. Thus, the condition of helicity conservation is most 
relevant in cases where SA = ^ x B is not exactly satisfied, 
i. e. resistive dissipation allows some motion of magnetic 
flux lines with respect to the fluid. Strictly speaking, such 
motion does not conserve either energy or helicity. However, 
helicity is more strongly dominated by large spatial scales 
than the magnetic energy, so small-scale resistive dissipa- 
tion may conse rve the form er to a better approximation 
than the latter (jFieldl . Il986f) . 

2.3. Toy model: a thin flux tube 

As a basis for later, educated guesses about the stability 
and evolution of MHD equilibria in stars, we examine the 
stability of a thin, azimuthal torus of cross section A lying 
in the equatorial plane of the star, at a distance r (^ VA) 
from the center, and containing a weak, roughly uniform 
azimuthal magnetic field B (^ ^/SnPt). For a general dis- 
cussion of th e properties of thin magnetic flux tubes, see 
Par ker] (jl979t ) and references therein. 

In order to be in equilibrium, the forces across the flux 
tube's cross section must balance, which requires the fluid 
pressure inside to be lower by SP = — _B^/(87r). This is 
achieved on the very short Alfven crossing time ^ A^^^/va, 
where va = B/^/Airph is the Alfven speed inside the flux 
tube. On a longer time ~ r/vA, the net forces on each sec- 
tion of the flux tube must also com e into balance. Its ten- 
sion, T = AB'^/{4tt) (jParkeii Il979[ ). causes a radial force 
per unit length fj- = —AB'^/iAirr) that tends to contract 
the flux tube. On the other hand, if the entropy and compo- 
sition of the matter inside and outside the flux tube are the 
same, the mass density inside will be lower than outside, 
Sp/pb — SP/{'-fpPi,) < 0, causing a radially outward buoy- 
ancy force per unit length, fg = —Spg, where g — iVi/^j is 
the gravitational acceleration. 

We take the point of view that the flux tube is initially 
placed at a radius r where the matter outside has the same 
composition and entropy as inside, and then allowed to dis- 
place to r -I- S^r, enforcing SP = —B'^/{87t) at each point, 
while the net force fnet = fx + fg controls the radial mo- 
tion. Using the notation of § [5J we note that 



Sp^ 



Pb 
IpPb 



<5P + ( J - 1 ) Pbg^r 



SO the net force per unit length can be written as 



fnet ~ 



A 
1^ 



1 - 



87r \7f, 



Pb9£.r 



(19) 



(20) 



where the pressure scale height H = Pb/\'^Pb\ = Pb/{pbg)- 
The term proportional to S,r accounts for the stratification 
of the fluid, and is manifestly stabilizing (force opposing 
displacement) if 7p > 7b and destabilizing (force reinforcing 
displacement) in the opposite case, while it vanishes for a 
barotropic fluid, for which 7^ — ^b- On the other hand, the 
term proportional to B^ is the force on the undisplaced flux 
tube, which will cause an inward displacement if r < 2jpH 
and an outward displacement in the opposite case, whereas 
r = 2jpH corresponds to an unstable equilibrium point. 



6 



Reisenegger: Magnetic equilibria in stars 



In a stably stratified fluid (7^ > 7;,), an equilibrium will 
be reached at the displacement 



1 - 



(21) 



Clearly, this equilibrium is stable with respect to radial, 
azimuthally symmetric displacements. However, it is intu- 
itive that the flux tube could contract towards the axis 
by moving away from the equatorial plane, roughly on a 
sphere of radius r. This motion would be driven by the 
tension, without being opposed by the buoyancy force. It 
could only be prevented by an additional, poloidal magnetic 
field, which can either be enclosed by the toroidal flux tube 
under consideration or be present in the form of a twist of 
the magnetic field in the tube. 

In all other cases (7p < 7f,), including that of a 
barotropic fluid (7^ = 7;,), there will be no equilibrium 
except at r = "^jpH, and the flux tube will either expand 
(if r > 2jpH) or contract (if r < 2jpH) indefinitely, at a 
speed determined by the fiuid drag force (Parker, 1974). 

This simple example suggests that, in the general case, 
the stratification of the fluid is likely to play an important 
role in determining the structure of magnetic equilibria, 
in the sense that there should be a much wider variety 
of possible equilibria in a stably stratified fluid than in a 
barotropic one. 

2.4. Axially symmetric equilibria 

The stable equilibria found by Br aithwa ite and 
coUaborators (Braith waite fc Spruitl . |2004 120061: 
iBraithwaite fc Nordl undl . |2006( ) ca n be described ide- 
ally as axially symmetric (but see IBraithwaite 2008i for 
highly asymmetric equilibria), involving two distinct 
regions: a thick torus fully contained in the star and 
containing a twisted toroidal-poloidal field combination, 
and the rest of space, containing a purely poloidal field 
that goes through the hole in the torus, and closing outside 
the star, in this way giving the external field an essentially 
dipolar appearance. It had long been speculated that such 
stable configurations might exist, but this has never been 
confirmed analytically (see IBraithwaite fc Nordlundl |2006| 
for a discussion of earlier work). 

For this reason, here we assume axial symmetry, allow- 
ing for both poloidal and toroidal field components. In this 
case, all the fiuid variables depend only on two of the three 
cylindrical coordinates, m and z. The most general, axially 
symmetric magnetic field can be decomposed into a toroidal 
component 

Bt = B{m,z)V(j) (22) 
and a poloidal component 

Bp = \/A{w,z) x\/(P (23) 

(|Chandrasekhar fc Prendergastl[l956D . This decomposition 
makes it explicit that the field depends only on two scalar 
functions, B and A, and explicitly satisfies the condition of 
zero divergence independ ently for both components. As in 
iReisenegger et all (|2007tl . we choose to write it in terms of 
the gradient of the azimuthal angle, V(/) = (f>/'cu, instead of 
the unit vector, 4>, in order to make easy use of the identity 



V X W(f> = 0. For reference, we also write the toroidal and 
poloidal components of the current. 



JT = X Bp = -^A*A 



4tt 



JP 



-^V xBt= -^-V6 X V0, 
An 47r 



(24) 
(25) 



where A* = ro^V • {zu ^V) is sometimes known as the 
"Grad-Shafranov operat or", although it appea rs to have 
been first introduced bv lLiist fc Schliiteil (jlpM). Eqs. ^ 
through (P5)) show that the magnetic field lines lie on the 
surfaces A = constant, while the current lines lie on sur- 
faces B = constant. If both A and B are taken to be zero on 
the symmetry axis, then 2TTA is the poloidal flux enclosed 
by a given surface A = constant, whereas cB/2 is the to- 
tal current en closed by the corresponding surface (see also 
lKulsrudll2005l § 4.9). 

In axial symmetry, the gradients in eqs. Q and ^ do 
not have an azimuthal component, and therefore eq. ([2]) 
requires ■ = 0, or equivalently 



jp X Bp = 0, 



(26) 



i. e. jp and Bp must be parallel to each other everywhere. 
(Note that jt and Bj- are always parallel.) In terms of the 
scalar functions defined above. 



yAx\/B=^ 0, 



(27) 



i. e. the surfaces A = constant and B = constant coincide, 
making it possible to write one of these functions in terms of 
the o ther, e. g. B = B{A) llChandrase khar fc Prendergastl 
Il956l ). When this condition is satisfied. 



1 



A* A- 



(28) 



with only two vector components (in the zu—z plane). Thus, 
for any choice of the functions A{zu, z) and B{A) whose 
magnitude is small enough to satisfy eq. (fTTj) . it should be 
possible to find independent scalar functions 6P and Ap in 
eq. (inj) t hat yield an eq uilibrium state. Thus, as already re- 
alized bv lMestell (|195Q) . any (weak) axially symmetric field 
satisfying eq. 1126]) corresponds to a magnetostatic equilib- 
rium in a stably stratified fluid. 

The possible equilibria are much more restricted in the 
barotropic case, in which the stabilizing Ap term in eq. ^ 
vanishes and the fluid force depends on a single scalar func- 
tion h = SP/pb. Using this, together with eq. (|28|) . in the 
force-balance equation one finds that V/i is parallel to 
\7A, so h = h{A), and eventua lly one obtains the popular 
Grad-Shafranov equation (e. g. lKulsrudll2005l § 4.9), 



A * ^ ^dB , ^ dh 



(29) 



which is often assumed to 
netic fields (Tomi mura fc Eri 
2006; HaskelletiS 



mag- 



characterize stellar 

uchi, 200,5 ; Yoshida et al 

Akgiin fc Wassermani . I200S 



.Kiuchi fc Kotakell2008( l. We emphasize that, in all the stars 
of interest here, the fluid is not barotropic, but stably strat- 
ified, with stabilizing buoyancy forces much stronger than 
the Lorentz forces, so the magnetic equilibria are not re- 
quired to satisfy eq. (|29|) , but only the condition contained 
in eqs. dHj) and (HTj). 



Reisenegger: Magnetic equilibria in stars 



7 



Of course, the existence of an equilibrium does not guar- 
antee its stability, which is clearly illustrated by the two 
simplest cases of purely toroidal and purely poloidal fields, 
for which there are equilibria, which however are always 
unstable to non-axisymmetric perturbations. For a purely 
toroidal field, flux rings can shift with respect to each other 
on spherical surfaces, in this way re ducing the total energy 
of the configuration (jTavleil . Il973f ). For a purely poloidal 
field, one can imagine cutting the star along a plane paral- 
lel to the symmetry axis and rotating one half with respect 
to the other, eliminating the dipole moment and reducing 
the energy of the external field, withou t changing the in- 
ternal one (|Flowers fc Rudermanl . I1977D . 

In the long-lived configurations found numerically by 
Braithwaite and collaborators, it is clear that the toroidal 
and poloidal field components might stabilize each other 
against both kinds of instabilities mentioned in the previ- 
ous paragraph. We can add here, based on the previous 
discussion, that the toroid of twisted field lines can be seen 
as a collection of nested, toroidal surfaces on which lie both 
the magnetic field lines and the current density lines (al- 
though their winding angles are generally different). As a 
consequence, the configuration has no toroidal Lorentz force 
component, although it generally does have poloidal com- 
ponents that are balanced by a pressure gradient and grav- 
ity 

We stress that Braithwaite's and his colleagues' simula- 
tions considered a single-fluid, stably stratified star. We can 
view the toroid, at least qualitatively, as a thick version of 
the thin flux tube discussed in ^ \2.'6\ It is impeded from con- 
tracting onto the axis by the presence of the poloidal flux 
going through it, as well as by the material between the 
torus and the axis. In the stably stratified star, the matter 
inside the toroid may have a slightly different entropy or 
composition than outside, cancelling its tendency to radial 
expansion due to buoyancy. However, if the star were not 
stably stratified (or this stratification could be overcome 
somehow; see § then a toroid sufficiently close to the 
surface would tend to rise and eventually move out of the 
star. If instead it were deep inside the star, it would natu- 
rally tend to contract due to tension, but be impeded from 
closing on its center by the poloidal flux. However, in this 
case, a small displacement of the whole configuration along 
the axis would cause a net buoyancy force that would tend 
to move it out of the star, along the axis. Although it is by 
no means clear whether this effect leads to an instability or 
instead is quenched by other effects, such as the progressive 
thickening of the toroidal ring or the material trapped by 
the poloidal field, we conjecture that stable equilibria occur 
only in stably stratified stars, and not in barotropic ones. 
If this conjecture were correct, it would make the usual 
search for barotropic (Grad-Shafranov) equilibria in fluid 
stars astrophysically meaningles^. 



3. Long-term field evolution 

From the previous discussion, it becomes natural to sug- 
gest that, as an a l ternative to the g e nerall y slow Ohmic 
diffusion ([Cowlind . 119451 : iBavm et al.l . Il969l ). the decay of 
the magnetic energy may be promoted by processes that 



^ They might, however, pl ay a role as stable "Hall equilib ria" 
in solid neutron star crusts (jLvutikov fc ReiseneggeJ . l2009f ). 



progressively erode the stable stratification of the stellar 
matter. 

For example, in an entropy-stratified, plane-parallel at- 
mosphere, a horizontal flux tube with a purely longitudi- 
nal magnetic field B can reach a mechanical equilibrium in 
which the interior entropy differs from the external one by 
Ss < 0, so as to compensate for the fluid pressure differ- 
ence induced by the magnetic field, 6P = —B'^/{8tt) < 0, 
yielding the same mass density inside as outside, Sp — 
{l/cl)[6P — {dP/ds)pSs] = 0, and thus zero net buoyancy. 
However, in this state, the temperature inside the flux tube 
is also lower than outside, so heat will strea m inwards, re- 
ducing I (5s I and making the fl ux tube rise (|Parkeil [l97l 
iMacGregor fc Cassinellil . |2003( ). This is the main alterna- 
tive to resistive diffusion in the case of entropy-stratified 
stars, i. e. upper main sequence stars and white dwarfs. 

Similarly, magnetic equilibria in a neutron star rely on 
a perturbation of the pr oton fraction, SY, which can b e 
reduced by two processes (jGoldreich fc Reiseneggeij . ll992[ ): 

1) direct and inverse beta decays, converting neutrons into 
charged particles (protons and electrons) and vice- versa, 
and 

2) ambipolar diffusion, i. e. diffusion of charged particles, 
pushed by Lorentz forces, with respect to neutral ones. 

In each case, if the magnetic structure was held in place, 
the imbalance {Ss or SY) would decay to zero on some char- 
acteristic diffusive or decay timescale tc- However, this de- 
cay corresponds to only a small fraction ~ i?^/ (SttP) = /3~^ 
of the absolute value of the relevant variable (s or Y), and 
therefore is compensated by a similarly small spatial dis- 
placement in the magnetic structure, ^ R/(3- A substan- 
tial change in the magnetic structure occurs only on the 
much longer time scale ~ /3ic- 

Since these processes and the corresponding tc differ 
substantially from one type of star to another, we now dis- 
cuss each type separately. 



3.1. Upper main sequence stars 

The Ohmic dissipation time for a magnetic field in a non- 
degenerate star is 



tn 



K 



3 X 10" yr 



/ T 



RqJ \10^K 



3/2 



(30) 



with Spitzer magnetic diffusivity rj — K/T'^/'^ and K ^ 
lO^'^ cm^ s~^ K'^/^. In order to obtain decay over the main- 
sequence lifetime of an A star, ~ 10^ yr, the characteris- 
tic length scale of the magnetic field configuration would 
have to be £ ^ 0.1 som ewhat smaller than found by 
iBraithwaite fc Spruit! (|2004[ ) . 

According to the discussion above, the time scale for 
decay of the field due to destabilization by heat exchange 
is ~ Ptc, where, in this case, tc is the heat diffusion time 
into a magnetic structure of characteristic scale I, related 
to the Kelvin-Helmholtz time, tKH, of the star (of radius 
R) by tc ^ (e/RftKH] therefore 



tB - P 



R 



> 1014 yr 



GAP 
RL 

£ 
R 



R 

Rq 



(31) 
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For realistic numbers, this time scale is comparable or some- 
what longer than the Ohmic time, thus not likely to be 
relevant for the star's magnetic evolution. 

3.2. White dwarfs 

In white dwarfs, the same processes are active as in main se- 
quence envelopes, although modified by the degenerate con- 
ditions. The Ohmic time scale is reduced (factor ^ 10^^) 
by the smaller length scale, and increased (factor ^ 10®) 
by the higher kinetic energy of the electrons (Fermi energy 
rather than thermal energy). Thus, again, the resistive de- 
cay of a large-scale field is too slo w to play a substant ial 
role in the evolution of these stars ([Wendell et al.l . [ToSTI) . 

Heat diffusion occurs chiefly through transport by 
the degenerate electrons, with conductivity k — 
STT^n^riek'^T/i'iZe^mfA), where h is Planck's constant, Z 
is the atomic number of the ions, e is the proton charge, 
TO* is the effective mass of the electrons (relativistic Fermi 
energy divided by , and A is the d imensionless "Coulomb 
logarithm" (see iPotekhin et al.l [199 S) , which we take ~ 1 
for the estimates that follow. Most of the heat is con- 
tained in the non-degenerate ions, whose number density 
is rii = Ue/Z, so the heat diffusion time[3 through a scale i 
is 



n.ke"^ 4 X lO'' yr / £ 



(32) 



Imposing the magnetic flux loss time ts ~ Ptc to be shorter 
than the coo hng t ime of the star, roughly given by Mestel's 
law ([Mestell . [T952[ ) as tcooi ^ 10^ yr/Tf 5, yields the condi- 
tion 



R 



< 



(33) 



Since /? > 10® and the temperature never drops below 
~ 10^ K, only very small-sc ale magnetic structures, ver y 
different from those found bv lBraithwaite fc SpruitI ([2004D . 
will be able to decay by this process. 

Unlike the case of neutron stars f? l3.3p . in known white 
dwarfs the thermal energy appears to be always larger than 
the magnetic energy, thus the eventual feedback of the mag- 
netic dissipation on the stellar cooling is negligible. 

3.3. Neutron stars 

Like white dwarfs, neutron stars are passively cooling ob- 
jects, in which the progressive decrease of the temperature 
makes the reaction rates and transport coefRcients (but not 
the spatial structure) change with time. In particular, with 
decreasing temperature, beta decay rates decrease dramat- 
ically, while collision rates also decrease and thus make par- 
ticle diffusion processes proceed more quickly. 

In the discussion below, we ignore the possibility of 
Cooper pairing of nucleons, which is expected to occur at 
least in some parts of the neutron star core and turns neu- 
trons into a superfluid and protons into a superconductor. 
This is likely to have a strong effect on the rates mentioned 



For white dwarfs, the diffusion time through scale ^ = i? in 
the degenerate interior is not the Kelvin-Helmholtz or cooling 
time, as the bottleneck for the latter is the conduction through 
the non-degenerate atmosphere. 



in the previous paragraph. However, this effect is difficult 
to quantify; therefore we rely on the better-known proper- 
ties of "normal" degenerate matter and leave it to future 
work to explore the Cooper-paired analog. 

3.3.1. Direct and inverse beta decays 

We first consider a hot neutron star (in the neutrino cool- 
ing regime, e. g. lYakovlev fc Pethicki [200l . in which the 
collision rates are so high as to effectively bind all particle 
species together, but weak interaction processes proceed at 
non-negligible rates. 

For illustration, let us again consider the toy model of 
a thin, horizontal magnetic fiux tube that is rising due to 
magnetic buoyancy through a degenerate gas of neutrons 
(n), protons (p), and electrons (e). Since the time scale to 
reach chemical equilibrium is much longer than any dy- 
namical times, the fiux tube can be considered to be in 
hydrostatic equilibrium, e. g. its internal mass density is 
equal to that of its surroundings, Sp = 0, and its internal 
fiuid pressure is reduced to compensate for the magnetic 
pressure, SP — —B'^/{8tt). These two conditions are only 
compatible if the fiuid inside the flux tube is not in chemical 
equilibrium, namely 



r] = fin - Pp- Me 



drj 
dP 



(34) 



where /Zi are the chemical potentials of the three parti- 
cle species. In order to change P without changing p, the 
composition, here parameterized by the proton or electron 
fraction, Y = rip/inn + Up) = n^jijin + rip), must change 
as well. 

For a simplifie d equation of state with no n-interacting 
fcrmions (sec, c. g. lShapiro fc Teukolskv^ll983^ ■ we show in 
Appendix [B] that {dP/dri)p « rig, so 



2.6 Bfg keV 



10^ 



^16 



K, 



(35) 



where we have assumed small perturbations, \r]\ fj,^, i. e. 
B ^ (87rne/ie)^/^ ^ 4x 10"'^^ G, easily compatible even with 
magnetar field strengths. (We took Ue ~ 2 x 10"^^ cm""^ for 
the numerical estimates.) 

In order for the fiux tube to move, \r]\ has to decrease 
by inverse beta decays, p -\- e n -\- Vf,., i. e. hy one of the 
same processes (direct or modified Urea) that control the 
cooling of the star. 

In the "subthermal" regime ([Haensel et al.l . [2003 ). |?7| < 
fcT, the available phase space for these reactions is deter- 
mined by the temperature, and the time scale for the decay 
of I77I is ~ 10 times shorter than that for the decrease of T 
([Reiseneggeil Il995[ ). On this time scale, rj would approach 
zero if the flux tube was held at its initial position. What 
happens is that, as Y is decreased by the beta decays, the 
pressure inside the flux tube increases, the tube expands 
and rises to find a new hydrostatic equilibrium in which it 
continues to be in a slight chemical imbalance as described 
by eq. ([55[) . This allows us to relate the logarithmic changes 
in the proton fraction inside the fiux tube and the temper- 
ature in the star as the latter cools and the fiux tube rises. 



dlnY 
d\nT 



10 



B^ 



{d-q/dY)p nfie 



(36) 
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Fig. 1. Magnetic field - temperature plane for a non- 
superfluid neutron star core. The dot-dashed horizontal 
lines show the initial temperature (just after core col- 
lapse), and the transition from neutrino-dominated (modi- 
fied Urea) to photon-dominated cooling. The dashed diago- 
nal line corresponds to the equality of magnetic and thermal 
energy. Above and to the left of the solid line, the star cools 
passively, on the time scales indicated in parenthesis along 
the vertical axis, without substantial magnetic field decay, 
so the evolution of the star is essentially a downward ver- 
tical line. Once the solid line is reached, magnetic dissipa- 
tion mechanisms become important and generate heat that 
stops the cooling. The subsequent evolution is expected to 
be roughly along this line, with temperature and magnetic 
field decreasing together, much more slowly than the pas- 
sive cooling. 

A substantial displacement of the flux tube corresponds to 
this quantity being > 1, i. e. it requires B > (n/Xe/5)^/^ ~ 
10^^ G, stronger than inferred magnetar fields and danger- 
ously close to violating the linear limit set above. (In addi- 
tion, it would require T 3> 10^^ K for consistency with the 
"subthermal" condition.) For weaker fields, the star cools 
too fast for the flux tube motion to keep up. 

On the other hand, in the strongly "suprathermal" 
regime, \r]\ > 5fcT, the induced inverse beta decays leave 
more thermal energy in the star than is emitted in the form 
of neutrinos, i. e. in the region in which this chemical im- 
balan ce is present, the Urea processes have a net heating 
efi^ect ([Fernandez &: Reiseneggeij . l2005D and might therefore 
be able to keep the st ar warm during a time long e nough 
for the field to decay ([Thompson fc Duncan[ . ll996D . This 
heating-cooling balance occurs at 



2 B 



16- 



(37) 



On this line (see Fig. 1), the thermal energy in the star, 

(38) 



is less than its magnetic energy, 
Eb ~ lO^" Ble erg, 
by a factor 



(39) 



(40) 



Et/Eb ~ 2 X 10-4 Tg, 

and therefore the cooling process of the star will be delayed 
by the inverse of this factor. 



3.3.2. Ambipolar diffusion 

At somewhat lower temperatures, collision rates are re- 
duced (due to the increased degeneracy and reduced num- 
ber of available quantum states) , allowing different particle 
species to drift with respect to each other. The Lorentz 
force only acts directly on the charged particles (protons, 
electrons, and perhaps others), pushing them through the 
neutrons. The magnetic flux is only frozen into the charged 
particle fluid, which moves through the neutral fluid as fast 
as the balance of Lorentz force and collisions allows. If the 
charged fluid contains only protons and electrons, it will be 
barotropic. If additional particle species are present, it will 
be stably stratified due to their density-dependent abun- 
dances; 

iGoldreich fc Reiseneggeij ([1992[ ) decompose the charged 
particle flux ricVc of ambipolar diffusion into two modes: 

1) Irrotational, with V • (ficVc) and V x (ricVc) — 0, 
which builds up pressure gradients in the charged par- 
ticle fluid, which need to be eliminated by weak inter- 
actions in order for the motion to proceed. 

2) Solenoidal, with V • (ricVc) — and V x {ricVc) ^ 0, cor- 
responding to an incompressible charged-particle flow, 
which does not cause pressure gradients and only needs 
to overcome the frictional force due to charged-neutron 
collisions. 

The solenoidal mode is analogous to the motion of a 
barotropic, incompressible fluid, which should be enough to 
overcome the constraints imposed by stable stratiflcation. 
In a n on-superfluid np e fluid, th is mod e proceeds on a time 
scale ([Goldreich fc Reiseneggeij Il992() 



^ambip ~ 3 X 10 



3 ^8^5 



causing a magnetic energy dissipation 



-Ef_ 



0.5 X 10^ 



erg s 



(41) 



(42) 



that can, at sufficiently low temperatures, balance the dom- 
inant cooling luminosity, be it neutrinos (here for the mod- 
ified Urea process). 



- 3 X lO^^rl erg s'S 



(43) 



or thermal electromagnetic radiation from the stellar sur- 
face, 



- 10^3 r|-2 erg s~\ 

The first will happen at Tg ^ B 
Ts ~ BiV' (see Fig. 1). 



(44) 

2 /5 

' and the second at 



3.3.3. Neutron star summary 

Strongly magnetized neutron stars appear to be subject 
to processes that can erode the stable stratification and 
therefore cause an MHD-stable field to decay on time 
scales shorter than their observable lifetime. These pro- 
cesses are weak decays, which are dominant at very high 
field strengths, and ambipolar diffusion, at somewhat lower 
field strengths. In both cases, these processes become im- 
portant only once the thermal energy in the star is sub- 
stantially less than the magnetic energy, and therefore the 
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latter acts as a large reservoir that keeps the star hot for 
much mo re than its c ooling time in the im-magnetized case 
(see also iThompson fc Duncan 1996; Pons et al. 2003). If 
the field decayed homologously, the star would evolve fol- 
lowing the line of heating-cooling balance in Fig. 1. In 
fact, the evolution is likely more complex, involving loss 
of st ability, followed by abrupt re-arrangements of the field 
(jBraith w aite fc Spruitl . 1200^ . but these effects should oc- 
cur roughly on the heating-cooling balance line. 

Of course, neutron stars also have a solid crust, whose 
elastic and yielding properties are still highly uncertain. At 
very high field strengths, the Lorentz forces will distort the 
crust, which might act essentially as a fluid. At lower field 
strengths, the crust might act as a valve, controlling the 
loss of magnetic flux. The relative importance of the decay 
mechanisms in the crust (Hall drift, crust cracking) and 
core is still unclear, depending on the uncertain properties 
of both. 



4. Conclusions 

This paper contains a general discussion of several physical 
issues related to the existence of large-scale, coherent mag- 
netic structures in upper main sequence stars, white dwarfs, 
and neutron stars. The main conclusions are the following: 



Appendix A: No force-free fields in stars 

Consider the following integral over a volume containing 
the star of interest: 

/ r • (j X B/c) dV^ <f nT.jdSj - [ TudV (A.l) 
Jv Js(v) Jv 

(e. g. iKulsrudI l2005l Chapter 4), where the Einstein sum- 
mation convention is being used, and the magnetic stress 
tensor is 



An Sir''- 



(A.2) 



The last term in eq. (jA.l|) is minus the total magnetic en- 
ergy within V, Ub = /y B^/{8n) dV > 0. The surface inte- 
gral, taking the surface to be a sphere of radius r, becomes 



— (b dVt r^B^cos{2(3), 
Stt / 



(A.3) 



where /3 is the local angle between B and f. Outside a 
star, B falls at least as fast as r~^, so this integral goes 
to zero as r ^ cx). Thus, the only way to have j x B = 
everywhere is to have Ub = 0, i. e. B = everywhere (ex- 
cept, perhaps, at a set of points of measure zero). This 
means that no magnetic stars can exist whose field is 
force-free everywhere. Th e "force-free" configurat i ons o f 
iPerez-Azorfn et all (|2006[ ) or iBroderick fc NaravanI (|2008f) 
do not violate this theorem, because they have current 
sheets with infinite Lorentz forces on the stellar surface. 



1) Magnetic forces in these objects are generally weak com- 
pared to pressure and gravity forces, and their matter is 
strongly stratified by entropy or composition gradients. 
This means that at least some components of the mag- 
netic forces can easily be balanced by other forces. Thus, 
there can be a wide variety of possible equilibria. These 
equilibria are not force-free; in fact, force-free equilibria 
are not possible in stars. 

2) If the magnetic structure is axially symmetric, the only 
constraint it has to satisfy to be balanced by pressure 
and gravity forces is that the azimuthal component of 
the Lorentz force must vanish. This means that there 
must be a set of magnetic surfaces of toroidal topology 
containing both the magnetic field lines and the current 
flow lines. Since the fluid is not barotropic, there is no 
need for the magnetic field to satisfy a Grad-Shafranov 
equation. 

3) It is difficult to give general criteria for stabil- 
ity. However, it is likely that, in a stably stratified 
star, poloidal and toroidal field components of similar 
strength could stabilize each other. In a barotropic fluid, 
it is possible that no stable equilibria exist, as the mag- 
netic field might rise buoyantly and be lost from the 
star. 

4) The long-term evolution of the magnetic field is likely to 
be governed by dissipative processes that erode the sta- 
ble stratification. Heat diffusion in main sequence stars 
and white dwarfs appears to be too slow to cause an 
observable effect over the life time of these stars. In 
strongly magnetized neutron stars, ambipolar diffusion 
and beta decays might be causing the magnetic energy 
release observed in magnetars. 



Appendix B: Thermodynamic properties of a 
degenerate npe fluid 

Taking the neutrons and protons to be nonrelativistic, the 
electrons extremely relativistic, and all highly degenerate, 
the total energy density and pressure are 

p = nmc^ + ^(n„/i„ 4- Upjlp) + ^Uefie, 



P 



n 



d{p/n) 



n 



'-{n„flr. 



n.. 



-nePe 



(B.l) 



Y 



where n„ = (1 — Y)n, Up — rie = Yn, /ie = 1ic{?ti:'^neY/^ , 
i^i — fJ-i ^ rnc^ = h'^[3n^ni)'^^'^ for i — n,p. The chemi- 
cal equilibrium state minimizes the energy per baryon with 
respect to Y at fixed n, so in equihbrium {dp/dY)n = 0. 
Thus, also evaluated at chemical equilibrium. 



dP_ 
9Y 

np., 
3 



2. ^ 1 

-{-Pn+ M + gMe 



The chemical imbalance rj — p^ 

— {l/n){dp/dY)n satisfies 



Pp P(^ — Pf^^ 



(B.2) 



drj 
W 



drj 
W 

Me 

3Y' 



2 Pn 

3 1 - y 



2pp 
3 Y 



1^ 
3 Y 



So, 



dP 
dp 



nY : 



(B.3) 



(B.4) 



Reisenegger: Magnetic equilibria in stars 



11 



Acknowledgements. The author thanks Stefano Bagnulo, Jon 
Braithwaite, Peter Goldreich, Swetlana Hubrig, Maxim Lyutikov, 
Friedrich Meyer, and Chris Thompson for many stimulating and in- 
formative conversations, Henk Spruit and an anonymous referee for 
insightful comments that improved the manuscript, and Cristobal 
Petrovich for preparing Fig. 1. This work was supported by Proyecto 
Regular FONDECYT 1060644 and Proyecto Basal PFB-06/2007. 



References 

Akgiin, T., & Wasserman, I. 2008, MNRAS, 383, 1551 
Arras, P., Cumming, P., & Thompson, C. 2004, ApJ, 608, L49 
Bagnulo, S., Landolfi, M., & Landi deglTnnocenti, M. 1999, A&A, 
343, 865 

Bagnulo, S., Landolfi, M., Mathys, G., & Landi deglTnnocenti, M. 

2000, A&A, 358, 929 
Baym, C, Pethick, C, & Pines, D. 1969, Nature, 224, 675 
Bernstein, L B., Frieman, E. A., Kruskal, M. D., & Kulsrud, R. M. 

1958, Proc. Roy. Soc, A244, 17 
Braithwaite, J. 2008, MNRAS, 386, 1947 
Braithwaite, J., & Spruit, H. 2004, Nature, 431, 819 
Braithwaite, J., & Nordlund, A. 2006, A&A, 450, 1077 
Braithwaite, J., & Spruit, H. 2006, A&A, 450, 1097 
Broderick, A. E., & Narayan, R. 2008, MNRAS, 383, 943 
Chandrasekhar, S., & Prendergast, K. H. 1956, Proc. Nat. Acad. Sci., 

42, 5 

Cowling, T. G. 1945, MNRAS, 105, 166 
Cutler, C. 2002, Phys. Rev. D, 66, 084025 

Donati, J.-F., Babel, J., Harries, T. J., Howarth, I. D., Petit, P., & 

Semel, M. 2002, MNRAS, 333, 55 
Donati, J.-F., Howarth, L D., Bouret, J.-C, Petit, P., Catala, C, & 

Landstreet, J. 2006, MNRAS, 365, L6 
Easson, I., & Pethick, C. J. 1977, Phys. Rev. D, 16, 275 
Fernandez, R., & Reisenegger, A. 2005, ApJ, 625, 291 
Ferrario, L., & Wickramasinghe, D. T. 2005a, MNRAS, 356, 615 
Ferrario, L., & Wickramasinghe, D. T. 2005b, MNRAS, 356, 1576 
Ferrario, L., & Wickramasinghe, D. T. 2006, MNRAS, 367, 1323 
Ferraro, V. C. A. 1937, MNRAS, 97, 458 

Field, G. 1986, in Magnetospheric phenomena in astrophysics, AIP 

Conf. Proc, 144, 324 (New York: American Institute of Physics) 
Flowers, E., & Ruderman, M. A. 1977, ApJ, 215, 302 
Goldreich, P., & Reisenegger, A. 1992, ApJ, 395, 250 
Haensel, P., Levenfish, K. P., & Yakovlev, D. G. 2002, A&A, 394, 213 
Haskell, B., Samuelsson, L., Glampedakis, K., & Andersson, N. 2008, 

MNRAS, 385, 531 
Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 
Hoyos, J., Reisenegger, A., & Valdivia, J. A. 2008, A&A, 487, 789 
Jones, P. B. 1988, MNRAS, 233, 875 
Kiuchi, K., & Kotake, K. 2008, MNRAS, 385, 1327 
KouveUotou, C, et al. 1998, Nature, 393, 235 

Kulsrud, R. M. 2005, Plasma Physics for Astrophysics (Princeton 

University Press) 
Lai, D. 1994, MNRAS, 270, 611 
Levin, Y. 2007, MNRAS, 377, 159 

Liist, R., & Schliiter, A. 1954, Zeitschr. f. Astrophys., 34, 263 
Lyutikov, M., & Reisenegger, A. 2009, in preparation 
MacGregor, K. B., & Cassinelli, J. P. 2003, ApJ, 586, 480 
Mastrano, A., & Melatos, A. 2008, MNRAS, 387, 1735 
Mathys, G., Hubrig, S., Landstreet, J. D., Lanz, T., & Manfroid, J. 

1997, A&AS, 123, 353 
Mestel, L. 1952, MNRAS, 112, 583 
Mestel, L. 1956, MNRAS, 116, 324 

Mestel, L. 1999, Stellar Magnetism (Oxford: Clarendon Press) 
Moss, D. 1984, MNRAS, 209, 607 
Moss, D. 1990, MNRAS, 244, 272 
Parker, E. N. 1974, Ap&SS, 31, 261 

Parker, E. N. 1979, Cosmical Magnetic Fields (Oxford University 
Press) 

Perez-Azorfn, J. F., Pons, J. A., Miralles, J. A., & Miniutti, G. 2006, 
A&A, 459, 175 

Pethick, C. J. 1992, in The Structure and Evolution of Neutron Stars, 
D. Pines, R. Tamagaki, & S. Tsuruta, eds., p. 115 

Petit, v.. Wade, G. A., Drissen, L., Montmerle, T., & Alecian, E. 
2008, MNRAS, 387, L23 

Pons, J. A., & Geppert, U. 2007, A&A, 470, 303 

Pons, J. A., Link, B., Miralles, J. A., & Geppert, U. 2007, 
Phys. Rev. Lett., 98, 071101 



Potekhin, A. Y., Baiko, D. A., Haensel, P., & Yakovlev, D. G. 1999, 

A&A, 346, 345 
Reisenegger, A. 1995, ApJ, 442, 749 
Reisenegger, A. 2001a, ApJ, 550, 860 

Reisenegger, A. 2001b, in Magnetic Fields across the Hertzsprung- 

Russell Diagram, ASP Conference Series, vol. 248, eds. G. Mathys, 

S. K. Solanki, & D. T. Wickramasinghe, p. 469 
Reisenegger, A. 2003, in Proceedings of the International Workshop on 

Strong Magnetic Fields and Neutron Stars, eds. C. Z. Vasconcelos, 

et al. (arXiv:astro-ph/0307133) 
Reisenegger, A. 2007, Astron. Nachrichten, 328, 1173 
Reisenegger, A. 2008, Rev. Mcx. Astron. Astrof., in press 

(arXiv:0802.2227 astro-ph]) 
Reisenegger, A., Prieto, J. P., Benguria, R., Lai, D., & Araya, P. A. 

2005, in Magnetic Fields in the Universe: From Laboratory and 

Stars to Primordial Structures, AIP Conference Proceedings, 784, 

263 

Reisenegger, A., Benguria, R., Prieto, J. P., Araya, P. A., & Lai, D. 

2007, A&A, 472, 233 
Reisenegger, A., & Goldreich, P. 1992, ApJ, 395, 240 
Reisenegger, A., & Thompson, C. 2009, in preparation 
Ruderman, M. 1972, ARA&A, 10, 427 
Schmidt, G., et al. 2003, ApJ, 595, 1101 

Shapiro, S. L., & Teukolsky, S. A. 1983, Black Holes, White Dwarfs, 

and Neutron Stars (New York: Wiley) 
Spruit, H. 2002, A&A, 381, 923 

Spruit, H. C. 2008, in 40 Years of Pulsars: Millisecond Pulsars, 
Magnetars, and More, AIP Conference Proceedings, 983, 391 
( arXiv: 0711.3650 astro-ph] ) 
Tassoul, M., & Tassoul, J.-L. 1983, ApJ, 267, 334 
Tayler, R. J. 1973, MNRAS, 161, 365 
Thompson, C, & Duncan, R. 1993, ApJ, 408, 194 
Thompson, C, & Duncan, R. C. 1996, ApJ, 473, 322 
Tomimura, Y., & Eriguchi, Y. 2005, MNRAS, 359, 1117 
Urpin, V. A., & Shalybkov, D. A. 1991, Sov. Phys. JETP, 73, 703 
Wasserman, I. 2003, MNRAS, 341, 1020 

Wendell, C. E., van Horn, H. M., & Sargent, D. 1987, ApJ, 313, 284 
Woltjer, L. 1958, Proc. Nat. Acad. Sci., 44, 489 
Woods, P. M. 1999, ApJ, 524, L55 

Yakovlev, D. G., & Pethick, C. J. 2004, ARA&A, 42, 169 
Yoshida, S., Yoshida, S., & Eriguchi, Y. 2006, ApJ, 651, 462 



